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Domain wall fermions (DWF) are a new fermion discretization scheme with greatly 
improved chiral symmetry. Our final goal is to study the nucleon spin struc- 
ture through lattice simulation using DWF. In this paper, we present our current 
progress on two topics toward this goal: 1) the mass spectrum of the nucleon 
excited states and 2) the iso- vector vector and axial charges, gy and <7a> of the 
nucleon. 



1. Introduction 

The RIKEN-BNL-Columbia-KEK QCD Collaboration has been pursuing 
the domain- wall fermion (DWF) methodEJ in lattice quantum chromodynamics 
(QCD). In DWF an extra fifth dimension is added to the lattice. By manip- 
ulating the domain-wall structure of the fermion mass in this fifth dimension, 
we control the number of light fermion species in the other four space-time 
dimensions. These light fermions possess exact chiral symmetry in the limit 
of an infinite fifth dimension. In particular: 1) fermion near-zero mode effects 
are well understoodcl, 2) explicit chiral symmetry breaking induced by a finite 
extra dimension is described by a single residual mass parameter, which is 
very small in the present calculation, in the low-enexgy effective lagrangianB, 
and 3) non-perturbative renormalization works welfi 

DWF is a promising new approach for treating fermions on the lattice. 
However, we need several tests of DWF in the baryon sector to reach our final 
goal of establishing the spin structure of the nucleon from first principles. 
Here we report our recent studies of the mass spectrum of the nucleon and 
its excited states and the nucleon matrix elements of the iso-vector vector 
and axial charges in quenched lattice QCD with DWF. Although most of the 
latter results are preliminary, the conclusive results in the former subject have 
been reported in Ref. [4]. 

2. Nucleon excited states 

First, we discuss the mass spectrum of the nucleon N and its excited states 
(the negative-parity nucleon N* and the positive-parity first excited nucleon 
N') by means of a systematic investigation utilizing two distinct interpolating 
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operators B^ 1 and B^ ■ For an explanation of those operators, see Ref. [4]. Our 
quenched DWF calculation was employed on lattice with size 16 3 x 32 x 16, 
gauge coupling (3 — 6/g 2 — 6.0, and domain wall height M5 — 1.8. Additional 
details of our simulation can be found in Ref. [4]. 

In Fig.l we show the low-lying nucleoli spectrum as a function of the 
quark mass, to/ in lattice units (a -1 w 1,9 GcV set from aM p =Q. 400(8) in 
the chiral limit). B^ gives the ground-state nucleon mass N (cross). The 
N* mass estimates (square and diamond) are extracted from both B± and 
B2 ■ The corresponding experimental values for N and N* are marked with 
lower and upper stars. Both N* mass estimates extracted from two distinct 
operators agree with each other. The large N-N* mass splitting is clearly 
evident. 

In contrast to the negative parity operators, we find that the mass esti- 
mates from a second operator, _B^~, are considerably larger than the ground 
state obtained from B± . This suggests that B^ has negligible overlap with the 
nucleon ground state and provides a signal for the positive-parity excited nu- 
cleon N'. To justify this possibility, we employ a sophisticated approach which 
utilizes the transfer matrix of a 2 x 2 correlation functional constructed from 
both Bf arid B% ■ The diagonalization of the transfer matrix yields the ex- 
cited stateQ. Fig. 2 shows a comparison of the fitted mass from ((B£ (t)B~£ (0))) 
(circle) and the estimated mass from the average effective mass given by the 
smaller eigenvalue of the transfer matrix (bullet). The cross symbol corre- 
sponds to the nucleon ground state mass evaluated from the larger eigenvalue 
of the transfer matrix, which is quite consistent with the fitted mass from 

((Btmtm 

Another important conclusion can be drawn from Fig.l and Fig. 2. In the 
heavy quark mass region, the ordering of the negative-parity nucleon (N*) 
and the positive-parity excited nucleon (N') is inverted relative to experi- 
ment. This remarkable result was originally reported in our early papeiu and 
subsequently confirmed in Ref. [6]. Further systematic calculation is required 
to determine whether this ordering switches to the observed ordering as one 
approaches the chiral limit. 

3. Nucleon matrix elements 

The nucleon (iso- vector) axial charge g A is a particularly interesting quan- 
tity. We know precisely the experimental value g A — 1.2670(35) from neutron 
beta decay. Why does g A deviate from unity in contrast to vector charge, 
g v =1? The simple explanation is given by the fact that the axial current is 
only partially conserved in the strong interaction while the vector current is 
exactly conserved. Thus, the calculation of g A is an especially relevant test of 
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the chiral properties of DWF in the baryon sector. In addition, calculation of 
g A is an important first step in studying polarized nucleon structure functions 
since g A — Au — Ad where (p, s\qf^^^qf\p, s) — 2s^Aqf with s 2 = —1 and 

We follow the standard practice!!! for the calculation of g v and g A . The 
three-point functions for the local vector current J? — qfj^q/ and axial cur- 
rent j{ = g/757iQ/ (i — 1, 2, 3) are defined by 

^%t')=Tr[P r J2(TB 1 (S,t)J^,t')B 1 (0,0))} (1) 

x,x' 

with r = V (vector) or A (axial) where P v = P + = (1 + ji)/2 and P A = 
P+lilb- For the axial current, the three-point function is averaged over i = 
1, 2, 3. The lattice estimates of vector and axial charges can be derived from 
the ratio between two- and three-point functions 



lattice _ Gp(t, t') - G r (t, t') 



G N {t) 



(2) 



where G N (t) = Tr[P + Y jS {TB 1 {x,t)B 1 {Q,Q))}. 

Recall that in general lattice operators 0\ a t and continuum operator O con 
are regularized in different schemes. The operators are related by a renormal- 
ization factor Zq: O con (fj,) = Zo(a[i)0\ a t{a)- This implies that the con- 
tinuum value of vector and axial charges are given by g r = Z r g 1 ^ ttlcc . In 
the case of conventional Wilson fermions, the renormalization factor Z A is 
usually estimated in perturbation theory (Z A differs from unity because of 
explicit symmetry breaking). For DWF, the conserved axial current receives 
no renormalization. This is not true for the lattice local current. An impor- 
tant advantage with DWF, however, is that—the lattice renormalizations, Z v 
and Z A1 of the local currents are the sameQ so that the ratio (3 A /<? v .) lattICC 
directly yields the continuum value g A . 

Our preliminary results are analyzed on 200 quenched gauge configura- 
tions at j3 = 6.0 on a 16 3 x 32 x 16 lattice with M 5 = 1.80. We choose a fixed 
separation in time of the nucleon interpolating operators, t = t sourc e — isink a 
t' < t, with currents inserted in between. We take t sourcc — 5 and £ s i n k = 21l 

In Fig. 3 we show the dependence of the vector renormalization, Z v = 
l/i?^ ttIC0 on the location of current insertions. A good plateau is observed 
so that Z v is certainly well behaved. The value 0.763(5) at m/ = 0.02 (ob- 
tained by averaging over time slices denoted by the solid line in Fig. 3) agrees 
well with Z A = 0.7555(3), which was obtained from a completely different 
calculation of meson two-point correlation functions based on the relation 
^conserved ( t )q l5 q(0)} = Z A (A 1 ™ 1 (t)q l5 q (0))i. 
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For the axial charge, g^ ttlce ; plateaus are seen for 10 < t < 16 in Fig. 4, 
so the charge ratios (g A /g v ) la,ttlce at each m/ are averaged in this time slice 
range. We find that there is a fairly strong dependence of (<? A /<?i,) lattICC on 
rrif as be shown in Fig. 5. A simple linear extrapolation to ai/ = yields 



simple linear ansatz may not describe the data in which case the result in the 
chiral limit may be even smaller. 

To compare to results using Wilson fermions, we plot g A versus the square 
of the ir-p mass ratio in Fig. 6. Our result is given by the bullets. We also 
include two heavier mass points (with large errors) from an earlier simulation 
using a larger separation between the source and the sink. The triangles 
and diamonds are from Wilson simulations at relatively strong! and weak 
coupling!! At first glance, the DWF and Wilson fermions seem to be in rough 
agreement except for the lightest point. However, our DWF results have a 
strong mass dependence. This may be a finite volume effect; the Wilson 
results at strong coupling were simulated on a lattice with a volume which is 
twice ours. 

A couple of comments on the mass dependence of our data are in order. 
First, it is interesting to note that our results appear consistent with the 
value 5/3 (marked as star) in the heavy quark limit, while the others seem 
inconsistent with this limit. Second, our results may also be consistent with 
vanishing in the chiral limit. This can be explained through the axial Ward- 
Takahashi identity which governs g A . If the PC AC relation cx mt is 
modified, for example by chiral logarithms, the nucleon matrix element of 
the pseudoscalar density does not develop a pole as m/ — > 0, and the r.h.s 
of the identity vanishes in the chiral limit. Thus, g A must also vanish as 
Tfif — > 0. Indeed, we already know that the PCAC relation for the pion mass 
is modified in the quenched approximation by two-effects: zero-modes of the 
Dirac operator and the quenched chiral logarithma. Further investigation of 
the Ward-Takahashi identity is under way. 

4. Conclusions 

We have explored several nucleon properties in quenched lattice QCD 
using domain wall fermions toward our final goal of studying the spin structure 
of the nucleon from first principles. 

Our quenched DWF calculation reproduces very well the largejnass split- 
ting between the nucleon iV(939) and its parity partner A*(1535)q. We have 
also calculated the mass of the first positive-parity excited state iV'(1440) by 
the diagonalization of a 2 x 2 matrix correlator and confirmed that it is heav- 
ier than the negative-parity excited state N* (1535)u. A remaining puzzle is 
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However, a 



whether or not a switching of N* and N' occurs close to the chiral limit. 

A preliminary calculation of iso-vector vector and axial charges shows 
that all the relevant three-point functions are well behaved. Z v determined 
from the nucleon matrix element of tha-vector current agrees closely with that 
from an NPR study x>f quark bilinearsu and a direct calculation using meson 
correlation functions^. This indicates that g v = 1 and Z v = Z A are mutually 
satisfied in our quenched DWF calculation. However, a linear extrapolation 
of g A to the chiral limit yields a value which is a factor of two smaller than 
the experimental value. We are currently investigating the Ward-Takahashi 
identity that governs g A to shed light on this behavior. We also plan to check 
related systematic effects arising from finite lattice volume and quenching (for 
example quenched chiral logarithms, zero modes, and the absence of the full 
pion cloud), especially in the lighter quark mass region. 
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Figure 1. N and N* (square and diamond) masses versus the quark mass my in lattice 
units. Note the large N-N* mass splitting which is within 10% (in the chiral limit) of the 
experimental value (bursts). 
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Figure 2. The mass of the positive-parity excited state (circles) is too high compared to 
the nucleon ground state (cross) to account for the observed splitting. 
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Figure 3. Dependence of vector renormalization, Zy = l/g^ ttlcc , on t', at rrif = 0.02. A 
good plateau is observed. 
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Figure 4. The lattice axial charge, g^ ttlco , at irif = 0.02. A good plateau is seen in the 
range 10 < t < 16. 
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Figure 5. Dependence of {g A / gv) l ^ ttlc ° on m/. 
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Figure 6. gA versus (m^/m p ) 2 . 
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